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Abstract 


Calculated response functions are presented which describe 
the pulse-height distributions produced by monoenergetic elec- 
trons incident on silicon detectors, It is assumed that the 
direction of incidence is perpendicular, and that the detectors 
are wide enough so that lateral leakage of energy is insignifi- 
cant, The validity of the calculations has been confirmed 
experimentally for incident energies between 0.25 MeV and 1.0 
MeV and detectors with thicknesses between 0.061 mm and 3.0 mm. 
The calculated response functions cover a wider range of condi- 
tions, including fourteen energies between 0.15 MeV and 90 MeV 
and ten detector thicknesses between 0.05 and 10.0 mm. At 
energies 2 1 MeV, response functions are given not only for bare 
detectors but also for an arrangement in which the detectors 
are shielded on the transmission side by another detector opera- 
ting in anti-coincidence, 

Key words: Electron; energy absorption; electron back- 

scattering; electron transmission; detector 


response function; silicon detectors, 
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Preface 


This report is an expanded version of a paper (published in Nuclear 
Instruments and Methods 6 5 181-193 (1969) ) on the response of silicon 
detectors to fast electrons, In that paper it was shown that transport 
calculations by the Monte Carlo method yield response functions that are 
in good agreement with experimental results, The new material in the 
present report consists of nineteen pages of tables listing energy 
deposition distributions D(T.,T) obtained in Monte Carlo calculations, 
These distributions, discussed in Section 2, relate the kinetic energy 
TS of the incident electron to the amount T of energy deposited in the 
detector, taking into account the transmission and reflection of energy 
carried by electrons and secondary bremsstrahlung, Tables of D(T,,T) 
are given for fourteen values of ses between 0.15 MeV and 5.0 MeV and ten 
detector thicknesses between 0.05 mm and 10.0 mm. These results can 
easily be converted to response functions by folding in the experimental- 
ly determined noise distribution (intrinsic resolution) according to the 


procedure given in the Appendix, 
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1. Introduction 


The use of silicon detectors for radiation measurements has become 
increasingly important in recent years. The interpretation of such 
measurements requires knowledge of the response of the detectors to 
various types of radiation, The purpose of this paper is to provide 
information about the response to electrons, 

Data on the response of silicon detectors to electrons have been 
reported by McKenzie and Ewan [1]*, West, Burns, Wang and Ziemba [2], 
Charoenkwan (3:I\, Rester and Rainwater [4], Huhn and Schneider eae 
Aitken, Emerson and Zulliger era Planskoy Bae and Waldschmidt and 
Wittig [8]. In the present investigation, extensive new data have been 
obtained experimentally as well as by a Monte Carlo calculation. The 
results provide information required for the application of silicon 
detectors to the counting of electrons and the measurement of electron 
spectra. They also confirm the accuracy of the Monte Carlo method used 
as a tool for predicting the penetration and diffusion of electrons in 
extended media. 

Experimental response functions and corresponding Monte Carlo 
results have been obtained for detectors with six thicknesses (0,061, 
0.105, 0.191, 0.530, 1.0 and 3,0 mm) at four electron energies (0.25, 
Os, O./5,-and 1,0 MeV) [9,10]. Monte Carlo response functions have also 
been calculated at fifteen energies between 0.15 and 5.0 MeV for 
detectors with ten different thicknesses between 0.05 and 10.0 mm, and 


are presented in Table 9 of the present report. 


ae 3 ne oss : 
Figures in brackets indicate the literature references at the end of 
this paper, 
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2e The Response Function 


The pulse-height distribution H(h), observed with a given silicon 
detector, is related to the actual electron spectrum S(T.) by the 


integral equation 


(ee) 


H(h) = J S(T RCT i) au a) 


The kernel R(T. yh) is called the response function’ of the detector, and 
R(T, h)dh is equal to the probability that an electron incident on the 
detector with kinetic energy Ty will give rise to a pulse with a height 
between h and htdh, 

It is useful to express the response function as the combination of 
an energy deposition distribution D(T.,T) and a noise distribution 


G(T,h) 5, by: writing 


T 
fe) 


R(T ,h) = | D(T,,T) G(T,h)ar. (2) 


In this expression D(T, ,T)dT is the probability that an electron incident 
with kinetic energy Ty will deposit an amount of energy between T and 
Td? in, the: detector, .'and G(T,h)dh is the probability that the deposition 


of energy T will give rise to a pulse with a height between h and htdh, 


2el Energy Deposition Distribution 


It is convenient to express this distribution as the sum of two 


terms 9 


DIT) P) DAT eae OSD, EL (3) 


A 
* The response function depends also on other parameters such as the 
detector thickness and the direction of the incident electron. In 
the present work, the direction of incidence is perpendicular to the 
detector surface. For data with oblique incidence, see Ref, 8. 


RQ 


where D(T, 57) is the probability density associated with the escape 


of some energy from the detector, P, is the probability of total absorp- 


A 
tion, and oe») is a delta function, The calculation of D(T,,T) was 


done by a Monte Carlo method which is described in Section 3, 
2.2 Noise Distribution 


The noise distribution, or intrinsic resolution, can be adequately 


represented by a Gaussian, 
G(T,h) =~ge— exp [- (hh )?/20? | (4) 
? 2-6 p Bais 


The most probable pulse height, pe is a linear function’ of T..° In this 

Section and in the Appendix we shall assume, without loss of generality, 
that the pulse-height scale is such that a is expressed in energy units 
and is equal to T.- The square of the full width at half maximum of this 


distribution is 
We = 8(£n2) 0°. (5) 


One can resolve We into the sum of three terms, Wi , We, and W3 , which 


take into account, respectively, (1) fluctuations in the number of hole- 
electron pairs produced, (2) fluctuations of losses occurring during 
charge collection, and (3) leakage current fluctuations and amplifier 
nOise.s 


In the present experiment, carried out at room temperature, We 


was found to have values ranging from 13 to 25 keV for different 


detectors, and to be independent of the incident electron energy. 


2 


¢ was WS which is the only 


Therefore, the predominant contribution to W 


term expected to be energy independent. The first term is relatively 


2 


small and is given by the equation 
Wy =8(fn2) FeT, (6) 


where € is the average energy required to peoduee a hole-electron pair, 
and F is the Fano [11] factor which accounts for the fluctuations of the 
number of pairs produced during the absorption of energy T. The experi- 
mental and theoretical information about these parameters has recently 
been reviewed by Deshpande [12] and Klein [13]. The latter author 
recommends the values € = 3.6 eV and F = 0.05 and thus obtains the 


1,0 MeV. The magnitude of the term WS can be 


value W, = 1 keV when T 
greatly reduced by increasing the bias voltage applied to the detector, 
Our data indicate that WS was small compared to Wes 


By combining Eqse (2-4), one obtains the response function in the 


form 
Is 2 2 
3 A -(h-T_)*/2o 
RUD ie Be Uh) acr caret ° ; (7) 
where 
= 
1 | 
h) === i le 
Re(T, yh) ae I dv.DA( 2.52) exp | (her) 2 /282 | | (8) 


The second term on the right-hand side of Eq. (7) is usually said to 


represent the absorption peak. R Toh) may be called the escape 


a 
component of the response function. The numerical evaluation of the 
integral in Eq. (8) is described in the Appendix for the case where a 


histogram representing DCT, »T) is available from a Monte Carlo 


calculation. 


3. Monte Carlo Calculation of the Energy Deposition Distribution 


3.1 Assumed Conditions 


A narrowly collimated monoenergetic beam of electrons was assumed 
to be incident perpendicularly ona silicon detector, For the purpose 
of the calculation, the detector was considered to be a plane-parallel 
target with lateral dimensions large enough so that the leakage of 
radiation through the side walls was insignificant. (This condition 
was also fulfilled in our experimental arrangement.) A one-dimensional 
treatment of the transport equation was therefore adequate, 

Most of our calculations were done for a single detector (case (i) 
im Fies 1), In addition, results were obtained for two special anti- 
coincidence arrangements (cases (ii) and (iii) in Fig. 1), by means of 
which the escape component of the response function can be substantially 
reduced, 


3.2 Sampling Procedures 


The following processes were taken into account: elastic scattering 
of electrons by atoms, inelastic scattering by atomic electrons and the 
production of knock-on electrons, the production of secondary bremsstrah- 
lung photons, and the scattering and absorption of these photons, A 
large number of collisions was sampled until the incident electron and 
the secondary electrons and photons were either absorbed in the detector 
or emerged from it, 

The sampling procedures have been described elsewhere rare and 
will be reviewed oniy briefly here. For the purpose of sampling, each 


electron track was divided into many short sections, chosen so that 


y) 


the electrons lost, on the average, a fixed fraction (~ 2%) of theta 
energy per section. dhe actual energy loss in each section, caused= by, 
many successive collisions with atomic electrons, was sampled from the 
Landau [16] energy-loss distribution (with the Blunck-Leisegang [17] 
correction to take into account binding effects), Each track section 
was further divided into four subsections which, in our Monte Carlo 
model, were assumed to be rectilinear, The electrons were allowed to 
change their direction at the end of each subsection, This change of 
direction was accomplished by sampling an angular multiple-scattering 
deflection from the Goudsmit-Saunderson | 18 | distribution, and com- 
pounding it with the original direction at the beginning of the sub- 
section. The Goudsmit-Saunderson multiple-scattering angular distribu- 
tion was evaluated with the Mott [19] elastic scattering cross section 
and a screening correction at small angles derived from the theory of 
Moliere [20]. Angular deflections associated with inelastic electron- 
electron collisions were taken into account approximately by replacing 
the factor Z* in the elastic scattering cross section by Z(Z+tl), 

where Z is the-atomic number of the medium, In each subsection, the 
following kinds of events were allowed to take place with the appropriate 
frequency: (a) the production of knock-on electrons, sampled according 
to the Mgller [21] cross section; (b) the production of bremsstrahlung 
photons, sampled according to a cross section based on the Bethe-Heitler 
theory with empirical modifications suggested in a review article by 
Koch and Motz [22]. The histories of the bremsstrahlung photons were in 
turn followed by conventional random sampling | 23]| through successive 


Compton scatterings until the photons either emerged from the detector 


6 


or were absorbed photo-electrically or in a pair-production event. The 
histories of all secondary Compton electrons, photo-electrons, or 


electron-positron pairs were followed in the same manner as those of the 


primary electrons. 
3.3 Utilization of Monte Carlo Histories 


For each source energy, the histories of 10* primary electrons, 
and the histories of all descendant secondary electrons and photons 
were sampled. This set of histories was analyzed to determine the 
energy deposition distribution for many detectors with different thick- 
NESSES» 

The sampled Monte Carlo histories were grouped into cascades, Each 
cascade consists of the history of a single primary electron and the 
histories of all of its descendant particles, Let f denote the 
fraction of the incident energy i that escapes through the transmission 
boundary of the detector, as the result of being carried away by primary 
or secondary particles, Similarly, let f. be the corresponding fraction 
of TS escaping through the reflection boundary, From each cascade, 
information was extracted pertinent to the three detector arrangements 
shown in Fig. 1: 

(aye “Ald cascades, regardless of the value of fo and Eps were 
allowed to contribute to the estimate of energy deposition for case ear 
The amount of energy left in the detector by the cascade was 
ice fe This energy was sorted into one of 50 energy bins covering 
the region from Ts to 0. 


(b) Cascades with fo = O were allowed to contribute to the estimate 


i 


of the energy deposition distribution for case (ii). The amount of 
energy deposited in this case was Ale eee 
(c) The relative frequency of cascades with ioe £,= 0 provided 


an estimate of the probability of total absorption, P 


A 


4e Experimental Method 


4.1 Electron Source 


The primary source of electrons was a d.c. accelerator which was 
stable to within + O.1 percent in high voltage. The absolute energy of 
the electrons was known to + 2 percent. In order to prevent pulse pile- 
up, a relatively small electron current was required, This requirement 


2 to 10m! 


was met by scattering the primary beam, which varied from 10 
amperes, from a thin gold foil. This foil was about 20 to 30 ug/em® 
thick, the gold having been vacuum-deposited on a 2- to 5-ug/em? thick 
film of collodion. As is shown in Fig. 2, the scattered electrons were 
carefully collimated at an angle of approximately 30 degrees with 

respect to the incident beam, and the electrons emerging from, theycol- 
limator were moving in a direction perpendicular to the face of the 
detector, The scattering chamber for this arrangement has been described 
elsewhere [9]. 

The scattered beam from the gold foil consisted mainly of elastic- 
ally-scattered electrons, but also contained some inelastically- 
scattered electrons of lower energy. The ratio of the inelastic to 
the elastic scattering cross section is an increasing function of the 
scattering angle. In order to minimize the relative amount of inelas- 


tically-scattered electrons in the beam, it was therefore desirable to 
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) 


make the scattering angle as small as possible. The angle of 30 degrees 
was the smallest that could conveniently be arranged with our experi- 
mental configuration. We have estimated that for primary energies 
between 0.25 and 1.0 MeV the amount of inelastically-scattered electrons 
in the beam at 30 degrees was of the order of one percent, The effect 


of this on the response function measurement is discussed in Section 5.l. 
4.2 Detectors 


Six silicon detectors which included surface-barrier, diffused, and 
lithium-drifted types were employed at room temperature. Some of the 
physical properties of the detectors are given in Table 1, For each 
detector the magnitude of the reverse bias applied was such that practi- 
cally the entire thickness was fully depleted, The dead layers at the 
surfaces were negligibly thin, except for the lithium-drifted detectors 
which had dead layers approximately 40 Um thick on their exit surfaces, 
According to the specifications supplied by the manufacturers, the 
detectors were uniform in thickness to l or 2 percent over the entire 
sensitive area. The area exposed to the electron beam was limited by 
collimation to a central region of 4 mm” , which was much smaller than 


the sensitive area, 
4.3 Amplification and Recording of Pulses 


The charge collected from the detector was fed into the input of a 
charge-sensitive preamplifier in which a pulse was formed, This pulse 
was further amplified and stored in 256 channels of a pulse-height 
analyzer, A rate of 400 to 800 pulses per second was maintained during 


each experimental run. The minimum height for pulse detection was 


a, 


established by a discrimination level (cut-off), hos which was set at 
the input of the analyzer. The values of ho adopted for the various 
detectors are given in Table l, 

The integral linearity of the entire electronic system was checked 
with a pulser and found to be within + 1 percent. The gain of the 
system was kept approximately the same for all detectors and source 
energies. By noting the position of the absorption peak of the response 
function at various source energies,” the linearity of the amplification 
system was confirmed, and the proportionality constant relating pulse 
height (channel number) and energy was established, With the gain 


setting employed, each channel had a width of 7.3 keV. 
4.4. Background Correction 


The background radiation incident on the detector was found to con- 
sist of bremsstrahlung photons produced in the beam trap and other parts 
of the apparatus, and a "dark current" of low-energy electrons originat- 
ing from components of the accelerator tube other than the filament, 

The background correction was obtained with two supplementary measure- 
ments: one with the scattering foil removed to determine the bremsstrah- 
lung background, and another with the accelerator filament cold and the 
scattering foil in place to determine the "dark current" background, 


The results showed that the bremsstrahlung contribution to the 


A check of the energy output of the d.c. accelerator was performed by 
matching the channel number vs energy relationship for the electrons 
from the accelerator to that obtained for conversion electrons from 
20°116(0.196 and 0.264 MeV), *87Cs5(0.626 and 0.656 MeV), and 
OU BIOLO7e and 1.044 MeV). This was done with a 3-mm thick detector, 
keeping the gain of the electronic system constant, 
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background was considerably larger than that of the "dark current", and 
increased with detector thickness. At primary energies of 0.5 MeV and 
lower, the entire background correction amounted to less than 1 percent. 
At higher energies, particularly for thick detectors (1 and 3 mm) the 
number of pulses in the background distribution amounted to 10 to 25 
percent of those in the principal distribution. The distribution of 
background pulses was concentrated in the vicinity of the cut-off value 
es decreased strongly with increasing h, and vanished for h approxi- 


mately equal to one-half the maximum pulse height. 
5. Results 


521 Comparison of Calculated and Measured Response Functions 


Response functions for various combinations of electron energy and 
detector thickness are plotted in Figs. 3a-d versus the relative pulse 
height, h/hoe The calculated results were obtained by folding the 
sampled energy deposition distribution into a noise distribution accord- 
ing to the procedure outlined in the Appendix. The width W,(see Fags 5) 
of the noise distribution, which was determined for each detector by 
fitting a Gaussian distribution to the experimental absorption peak, 
was found to be independent of the electron energy for each detector 
and is given in Table 1. Each experimental response function shown is 
based on a sample of Be tOe counts. The electron current incident on 


the silicon detectors was not monitored, Nevertheless, it was possible 


1 The quantity h_ is the most probable pulse height which is associated 
with the deposition of energy Ty in the detector, When pulse heights 
are expressed in energy units, then ho = Toe 


Lk 


to normalize the experimental response functions absolutely, and 
independently of the Monte Carlo calculations, by requiring that the 
total integral ) R(T. yh)dh be equal to unity. Actually, the pulse-height 
distributions nere obtained only above the cut-off value hoe However, 

in most cases the response function rapidly became small with decreasing 


h, so that the additional contribution to the integral from pulses 


’ 
be low he was negligible. In a few exceptional cases, namely, for the 
O.,061-mm detector at all re and for the 0.105-mm detector at 
Dee 0.50 MeV, the normalization was performed by setting the area under 
© 

the experimental response function down to the cut-off, ii R(T. ,h)dh, 
equal to the detection efficiency Na obtained in the Non cemen ere calcu- 
lations (see Eqe 9 of Section BN 

The agreement between calculated and experimental response functions 
is generally good. For certain cases with energies of 0.25 and 0.50 
MeV in Figs. 3a-d the experimental peak of the escape component occurs 
at a somewhat smaller relative pulse height than the corresponding cal- 
culated peak. It appears that in these cases the most probable energy 
loss of the electrons transmitted through the detector is somewhat 


smaller than predicted, We have been unable to explain the magnitude 


and pattern of these discrepancies in terms of experimental errors or 


1 In some of the experimental response functions, a small peak can be 
seen which is produced by inelastically scattered electrons from the 
gold foil (see Section 4,1). This peak appears at a relative pulse 
height h/h, ~ 0.6 as predicted by the Méller formula for electron- 


fo) 
electron scattering at 30 .« 
12 


statistical errors in the Monte Carlo calculations. The Monte Carlo 
model does not take into account "channeling effects", presumably caused 
by diffraction, of the type found by various authors [26-29] in electron 
penetration experiments in crystalline media. It is possible that this 


omission is responsible for the discrepancies. 
5.2 Scaling 


The extent to which a detector reflects or transmits electrons 
depends primarily on the ratio of the detector thickness, z, to the 
mean range,” roy of the incident electrons. A list of range values for 
silicon is given in Table 2. When the ration 2/ is small, the absorp- 
tion peak is either absent or quite small, and there is a large escape 
“component associated with the transmission of electrons. As 2/T 
increases, the absorption peak becomes more prominent, and the escape 


component is associated mainly with the reflection of electrons. The 


TThe random numbers needed for the Monte Carlo calculations were 
obtained with a multiplicative congruential generator, It has recently 
been shown by Marsaglia [25] that such generators have a defect which 
makes them unsuitable for certain Monte Carlo problems, The statistical 
analysis of different runs made with our computer program, and the 
comparison of results obtained with different congruential generators 
used singly or in combination, produced no evidence indicating a lack 

of randomness which would tend to bias the results, In particular, it 
appears very unlikely that the observed discrepancies between calcula- 
tion and experiment are due to a deficiency of the random numbers used. 


“The mean range is calculated from the expression 


ih 
os ae at/(- =) ; 


dE. 
where ax is the mean energy loss per unit pathlength (stopping power) 


due to collision losses and bremsstrahlung, 
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shape of the response function, for a given value of 2/T is very insen- 
sitive to the value of the incident energy Tos This can be seen by 
inspection of Figs. 3a-d and is brought out even more clearly by two 


experimental examples in Fig. 4. 
503 Backscattering, Transmission and Absorption Coefficients 


The response function calculations also give values of various 
reflection and transmission coefficients for plane-parallel silicon tar- 
gets. These coefficients may be defined for primary (incident) electrons, 
secondary electrons (delta rays, photoelectrons, Compton and pair elec- 
£rons), and bremsstrahlung. They are in each case equal to the average 
number of particles per incident electron that emerge from the entrance 
(backscatter) surface or the exit (transmission) surface. Coefficients 
for l-mm and 3-mm thick detectors are shown in Table 3. It can be seen 
that at high energies, because of the escape of secondary particles, the 
probability of total absorption of an entire electron-photon cascade is 
significantly smaller than the probability of absorbing the primary 
electron. 


The probability of total absorption, P for a large number of 


A? 
different detector thicknesses and incident energies is given in Table 4. 
The estimated uncertainty in these results is + one unit in the last 
figure, The statistical error is given by ye (By oe , where 105 cs 
the size of the Monte Carlo sample, The systematic errors include the 
effect of round-off, approximations in the Monte Carlo model, and un- 


certainties in the cross sections, and are believed to be of the same 


order of magnitude of the statistical error, 
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Values of the fraction of events without total energy absorption, 
CEP) reported by various authors are compared in Table 5. The 
present results (experimental and calculated) and Planskoy's results |'7| 
are in good agreement, They show a moderate rise with increasing qT, 
which is to be expected because of the increased escape of secondary 
particles, By contrast, the reflection coefficient shows a decrease 
with increasing Toe The values given by Charoenkwan bee Rester and 
Rainwater Wale and Huhn and Schneider [5] are higher than the present 
results at all energies, and those of Waldschmidt and Wittig [8] at some 
energies, The experimental problems involved in the determination of 
(lea) have been carefully investigated by Planskoy who found that, 
depending on the detector and collimator geometry,spuriously high values 
can easily be obtained in the following circumstances: (a) when the 
incident electrons can lose some of their energy before reaching the 
detector, as the result of inelastic: scattering in the collimator, and 
(b) when the electrons can escape with some energy through the side 
walls of the detector, In the present experiment, precautions were 
taken to make these effects negligibly small,” 

Table 5 also gives comparisons of the fraction of electrons back- 
scattered, The calculated values are in good agreement with the exper- 


imental result of Rester and Rainwater [4] at l MeV, but are one-third 


1 

It was necessary, however, to lower the measured values of (1-P,) 
by’ ~ 0,005,:.in order to.discount the presence of inelastically 
scattered electrons in the incident beam from the gold scattering 
foil. 
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smaller than the coefficients calculated by Waldschmidt and Wittig [8] 

at 0.3 and 0.6 MeV. The predictions of our Monte Carlo program have also 
been compared with many other related experiments, For example, the 
fraction of electrons backscattered from aluminum at 0.6 MeV has been 
calculated ‘to be 0.105. +°0,003, which is in good agreement, with the jex- 


perimental values of Glazunov and Guglya [30] (0.094), Cohen and Koral 


[any ¢ +0 .00'/ +0 .002 


0.104 0.002)? and Miller and Hendricks [ 32] (0.109 6°91) 


54 Detection Efficiency 


In’practice, only those pulses are reeordéd that are greater than 
some chosen threshold value, hay in order to reduce the pulses produced 
by noise and background radiation. The detection efficiency, Na may 
then differ significantly from unity, and can be calculated as an in- 


tegral over the response function, i.e., 


It is not practical to make a complete tabulation of the efficiency 
because it depends on a large number of parameters (e.ge., incident 
energy, detector thickness, threshold, width of noise distribution). fin 
Tables 6a-b, sample results are given which indicate the dependence on 


the various parameters, 
525 Anticoincidence Arrangements 


The escape component of the response function can be reduced through 
the use of an arrangement in which a detector operating in anticoinci- 


dence is placed on the transmission side of the principal detector 
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(case (ii) in Fig. 1). The extent of the reduction is demonstrated in 
Fig. 5 for incident energies between 2 and 5 MeV and a 3-mm detector, 
The efficiency 1h of this arrangement, i.e., the fraction of events in 
which an incident electron gives rise to a recorded pulse, is shown in 
Table 7 for various incident energies and detector thicknesses. In case 
(aaa) with the silicon detector surrounded on both sides with detectors 
operating in anticoincidence, the efficiency iN is, of course, equal to 
the probability of complete absorption, Pas It is also of interest to 
consider the ratio B/N, which represents the fraction of the area of 
the response function that is under the absorption peak. As can be 

seen in Table 8, the values of this ratio, for incident energies greater 


than 1 MeV, are considerably greater in case (ii) than in case (i). 
5.6 Tables of Energy Deposition Distributions 


The Monte Carlo calculations provide estimates of the energy depos- 
ition distribution in the form of histograms, The escape component of 
the distribution in histogram form is given by 


uf 


_T re [ n+1) 


n+l n T Put 
n 


psn) (T.) = (T T TAT, n=0,1,2e0. (10) 


where TL and Tel are the boundaries of the nt+l'st energy interval, The 
histograms are given in Table 9 for all of the incident energies Ty 

and detector thicknesses z listed in Table 4. For I> 1 MeV, results 
are also given for the anticoincidence arrangement (4a); with an 


auxiliary detector operating in anticoincidence on the transmission side, 
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Appendix 


Convolution of the Energy Deposition and Noise Distributions 


Given the energy deposition distribution in histogram form as 
listed in Table 9, and knowing the experimental value of the standard 
deviation 6 of the noise distribution, one can easily determine the 
corresponding response function. The histogram of the escape component 


of the response function may be defined, in analogy to Eq. (10), as 


ai") (Tr) = (nag) Ra(T,,h)dh. (11) 


The energy deposition and response function histograms are then related 


by a noise matrix G, i.e. 


(11) a om) 
=o) aD emer. (12) 
m 
where 
men Oe Wo pial gy fmt eax (eho Teen (13) 
mn n+l on i he Jot o 


The evaluation of the double integral in Eq. (13) gives the result that 


G = iG eee if r/ Anim) F( “nn is 
ia Cae cee un cy 2 Of. 
h =) h ae ae 
+1 m1 ("ttl m\ ; 
(ele, (el ee 
\ 3/2 bf2 / 


2 
where F(x) = x erf (x) + = iis 


x 2 


and erf (x) == i e © dt is the error function. 
re) 
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Table 1. Characteristics of the Silicon Detectors 


t 

| Detector eee oe Low-energy 
| Thickness ; Area Cutoff, h 
| (nm) | (mm) 


Surface barrier 


Ww 


Lithium drifted 


| 
| 
| 
Diffused | 
| 
| 
| 


* 
At room temperature 


on 


Table 2. Mean range .of ,electrons in silicon* 


ee cee oy 


eee _ : 
(MeV) (mm) | 
5.00 12s | 
4.00 1029 
3.00 Thales! 

2.00 5.02 

1.50 3.67 

1.00 2.29 

Once 1.74 

0.60 120 | 
0.50 0.936 

0.40 0.687 

250 Oei5 

0.25 Ons ut 

0.20 3 0.2h2 : 
0.15 3 0.152 | 


* 
Calculated by multiplying the values (in 
g/ cnt ) in Ref. 24 by 10/p, where 


9 = 2.55 g/em as the density, of silicon. 
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Table 6. Detection Efficiency, Nap as a Function of Various Parameters 


(a) Dependence on Electron Energy To) Cut-Off nergy ho» and Detector Thickness z 


Ng 
T (MeV) h_ (keV) 

z (mm) 0.5 TO) 2.0 0.5 a 20 
h.0 100 0.98 0.98 0.99 0.95 0.98 
200 0.17 0.97 0.99 0.19 0.96 
300 0.05 0.68 0.98 0.05 0.67 
Tele) 0.02 0.20 0.98 0.02 0.22 
O13 0.5 
2.0 0.9) 0.99 
0.36 0.96 
0.10 0.62 
0.05 0.23 
Onl 0.2 
1.0 0.83 0.99 
0.29 0.92 
0.08 0.60 
0.03 0.27 
0.05 Owal 
0.5 0.80 0.98 
O.4y 0.91 
0.17 On 7a: 
0.07 0.5 
0.05 0.1 
0.25 0.93 1.00 
OR TAL 0.99 
0.42 0.97 
0.23 0.91 
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Table 9. Energy deposition distributions calculated by the Monte 
Carlo method. The results are given in histogram form, The first 
two columns on the left give the limits a ye energy intervals, 
and the other columns give the quantity D\"/(T_) defined by Eq. (10). 
The numbers in the bottom row, labelled "' OTAL ABSORPTION", are 
identical with the values of P, given in Table 4. The numbers fol- 
lowing the symbols (+) or (-) indicate powers of ten, 
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Fig, 1, Arrangements for silicon detectors considered in calculations: 
Case (i): Single detector; Case (ii): Detector 2 in anticoincidence 


with detector 1; Case (iii): Detectors 2 and 3 in anticoincidence with 
detector l, 
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Fig. 2. Experimental arrangement of scattering foil, collimator, and 
silicon detector, 
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Fig. 3a. Response functions of silicon detectors for various thick- 
nesses Z To = 0.25 MeV. The points are experimental and the solid 
lines are calculated. 
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Fig. 3b. Response functions of silicon detectors for various thick- 
nesses Ze To = 0.50 MeV. The points are experimental and the solid 
lines are calculated. 
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Fig. 3ce Response functions of silicon detectors for various thick- 
nesses Z. To = 0.75 MeV. The points are experimental and the solid 
lines are calculated. 
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Fig. 3d. Response functions of silicon detectors for various thick- 
nesses Zz. Tp = 1.00 MeV. The points are experimental and the solid 
lines are calculated, 
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Fig. 5. Effect of an anticoincidence arrangement on calculated response 
functions, Case (i): Single detector; Case (ii): Single detector 
with an auxiliary detector operating in anticoincidence, 
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